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Abstract 

At nanohertz frequencies gravitational waves (GWs) cause variations in time-of-arrival of pul- 
sar signals potentially measurable via precision timing observations. Here we compute very-low- 
frequency GW sensitivity constrained by instrumental, propagation, and other noises fundamen- 
tally limiting pulsar timing observations. Reaching expected GW signal strengths will require es- 
timation and removal of ~99% of time-of-arrival fluctuations caused by typical interstellar plasma 
turbulence and a reduction of white rms timing noise to ~100 nsec or less. If these were achieved, 
single-pulsar signal-to-noise ratio (SNR) = 1 sensitivity is then limited by the best current ter- 
restrial time standards at h^-ms ~2 xlO^"*^^ cycle/year)] ^^/^ for / < 3 x 10^^ Hz, where f is 
Fourier frequency and a bandwidth of lcycle/(10 years) is assumed. This sensitivity envelope may 
be optimistic in that it assumes negligible intrinsic pulsar rotational noise, perfect time transfer 
from time standard to observatory, and stable pulse profiles. Nonetheless it can be compared to 
predicted signal levels for a broadband astrophysical GW background from supermassive black hole 
binaries. Such a background is comparable to timekeeping-noise only for frequencies lower than 
about 1 cycle/(10 years), indicating that reliable detections will require substantial improvements 
in signal-to-noise ratio through pulsar array signal processing. 

PACS numbers: 98.80.-k, 95.36. +x, 95.30.Sf 
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The pulsar timing GW detector uses the earth and a distant pulsar as electromagnetically- 
tracked separated test masses. The pulsar emission serves as a clock, in the idealized case 
producing perfectly periodic radio pulses transmitted to the earth. These are timed by 
cross correlation of the received pulses against a template of the pulsed waveform. Time-of- 
arrival residuals, R(t), are produced after correcting for known effects. Signals and noises 
enter the observed time series via transfer functions [Ij. A GW of characteristic strain 
amplitude h incident on the earth-pulsar system produces variations of order h in the time 
series of relative dimensionless frequency fluctuations, y(t), of the pulsar signal [2] y(t) = 
((1 - fi)/2)[^{t - Ti(l + /i)) - ^(t)]. Here /i = k • n, k is a unit vector parallel to the GW 
propagation direction, n is a unit vector from the earth to the pulsar, h(t) is the GW strain 
tensor, \E'(t) = (n ■ h ■ n)/(l — /i^), and Ti is the one-way light travel time between the 
pulsar and earth. The fractional frequency time series is the derivative p] of the observed 
time-of-arrival residuals, R(t): y(t) = dR(t)/dt. The pulsar timing technique has been used 
to bound GW signal strengths, e.g. [IJ-E]. 

To assess instrumental and other noises currently and fundamentally limiting detections, 
we compute here the sensitivity of pulsar GW observations. Sensitivity is conventionally ex- 
pressed as the sky- and polarization-averaged sinusoidal signal strength necessary to achieve 
a given signal-to-noise ratio in a given bandwidth, e. g. [21 [10] ■ Exphcitly, we compute 
the signal strength required to produce SNR = 1 in bandwidth B: [S'j,„(/)i?]^/^/(rms signal 
response), where Syn{f) is the spectrum of noise and the rms signal response in general also 
depends on Fourier frequency. 

The GW signal response depends on the angle of arrival of a wave relative to the earth- 
pulsar line. For GWs from a specific direction the above formula for y(t) can be used 
directly j^. We are interested here in signal response averaged over the sky. To get the 
rms signal response as a function of Fourier frequency, the Fourier transform squared of 
the GW signal response, above, is averaged over the sky and polarization states to obtain 
Sy{f)/Sh{f) = 1/3 - l/iSn^pTi'^) + sin(47r/Ti)/(32/37r3Ti2), where Sy is the spectrum 
of fractional frequency fluctuations, Sh is the spectrum of GW strengths, and f is Fourier 
frequency. In the practical case Ti is hundreds of years or longer and the duration of 
pulsar timing observations is ~decades, so the second and third terms are negligible for 
/ > l/(duration of observations): Sy{f)/ Sh{f) ~ 1/3, implying that rms signal response is 
constant 0.58) over the accessible frequency band. 

2 



To compute sensitivity we also need spectra of the noises. Important noise sources in- 
clude finite signal-to-noise ratio in the raw observations, instability of the local clock against 
which pulsars are timed (and errors in time transfer if the clock is not located at the ob- 
servatory), uncertainties in solar system ephemerides (used to correct arrival times at the 
earth to the barycenter of the solar system), pulsar position uncertainty, intrinsic pulsar 
rotational stability, stability and accuracy of the average pulse templates used to measure 
R(t), dispersion measure (DM) variability in the interstellar and interplanetary plasmas, 
tropospheric scintillation (the wet and dry components of the troposphere cause delay vari- 
ability), antenna mechanical noise (stability of the phase center of the antenna tracking the 
pulsar), and station location errors (changes in antenna location due to atmospheric and 
tidal loading of the crust). Figure ([T]) shows the GW sensitivity for several of these noise 
sources individually, as discussed briefly below, with sensitivity computed for bandwidth B 
= 1 cycle/10 years. 

The green curve is the station- location- noise limit; Syn was computed from 30 years of 
absolute value of vector ground displacement using data ^Jj taken near the NASA/JPL 
Goldstone CA tracking complex. The derivative theorem for Fourier transforms [12j was 
used to convert the spectrum of displacement to the spectrum of velocity and hence the 
spectrum of y = Av/c. The black curve is derived from the power spectrum of hourly zenith 
dry tropospheric pressure fluctuations [13] provided by the National Climate Data Center 
and taken at the NASA/JPL deep space tracking complex near Madrid, Spain. Pressure was 
converted to zenith path variation using (path variation in centimeters) = 0.022768*(surface 
pressure in millibars), ignoring a factor close to unity which depends on latitude and height. 
The blue curve is from the power spectrum of zenith wet tropospheric path delay, computed 
from 1.5 years of data [14j taken at the NASA/JPL Goldstone, California tracking com- 
plex. Tropospheric path variation spectra were similarly converted to spectra of y using the 
derivative theorem. The light blue curve is the measurement (for f > 0.0001 Hz) and upper 
limit (for 10~^ — 10~^ Hz) for antenna mechanical noise fluctuations observed with a 34m 
tracking antenna at Goldstone [13]; smaller stiffer antennas give lower antenna mechanical 
noise [16] [17] . The solid black lines are for white timing noise with rms amplitudes of 100 
nsec, 1 nsec, and 1 ps in a Fourier band ± 0.5 cycles/day (i.e. one sample per day; current 
observations are more typically 1 sample per 2 weeks, which would result in curves \/l4 
higher). The 100 nsec level is the current timing goal of leading timing array experiments; 
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three pulsars are being timed to this level [9]. One picosecond is the absolute best possible 
timing accuracy one can achieve using millisecond pulsars. Since the pulsar signal itself is 
an amplitude modulated noise process, it can be said to have self-noise. In the absence of all 
other sources of noise, including timing and antenna noise, the pulsar signal self-noise would 
still be present. Assuming that the narrowest possible average pulse profile is 10 /isec, the 
pulsar signal bandwidth is 1 GHz, and that one can observe the pulsar for 12 hours at a 
time, the self-noise yields an rms timing accuracy of about 1 ps. The dotted curve shows 
approximate limits due to the uncertainties in the masses of the planets P^] - EU] . Mercury 
through Jupiter, affecting knowledge of the solar system barycenter. 

The dashed lines are noise from representative interstellar medium dispersion measure 
(DM) variation [2T1 [22] (assumed Kolmogorov spectrum with = 0.001 m~^°/^, propa- 
gation distance z = 1 kpc, radio frequency = 1 GHz, transverse velocity v = 100 km/sec), 
SynW) = 7r-^^^2-^/h^/^X'^c-^zCn^re^[T[4/3]/r[ll/6]]f-^/^ with the indicated levels of cal- 
ibration, i.e. 99% calibration means only 1% of the DM fluctuation rms noise remains 
in the measurement. (Some nearby pulsars have smaller integrated ISM turbulence levels 
and would require smaller percentage corrections to reach the indicated line in Figure ([T])). 
The effect of solar wind plasma turbulence is non-negligible but its dispersive character 
should allow it to be calibrated in addition to the interstellar plasma. 

Pulsars are timed against terrestrial clocks. Recent stability measurements of linear ion 
trap time standards [21] give Syn{f) ~ 4 x 10~^^{f /lHz)~^ , measured in the approximate 
band 10^'' — 10^^ Hz. We assume this spectrum continues to be valid to lower frequencies; 
in a 1 cycle/ (10 year) bandwidth this noise gives GW sensitivity shown as the dot-dash line 
in Figure The dotted curve in the lower right, for comparison, is the low- frequency 
segment of the Laser Interferometer Space Antenna (LISA) missions predicted sensitivity 
curve ([25], S-cr- in-one-year- integration). 

Not included in Figure ([T]) are intrinsic pulsar rotational instability noise [2S1 [2Z] (variable 
by pulsar and substantial for some), errors due to time transfer from the frequency standard 
to the observatory p8], pulsar position uncertainty ^J, errors in pulse templates, and errors 
due to radio frequency dependence and temporal instability of pulse profiles. So the upper 
envelope sensitivity in Figure ([T]) is in this sense optimistic. 

One application of our sensitivity analysis is the detectability of an astrophysical GW 
background from incoherently radiating supermassive black hole binaries. Such a back- 

4 



ground is predicted to produce signal strengths [Uj I2UH5^ in the range h{f) = [fSh{f)Y^'^ ~ 
(1— to— 10) X 10~^^[f/ (1 cycle/year) ]~^/^. For comparison with our SNR = 1 in a fixed B = 1 
cycle/(10 years) bandwidth, we convert to Kmsif) = [BSh{f)Y^'^ ~ (3-to-30) x IQ-^^ [f/(l 
cycle/year)] These GW strengths are comparable to the SNR = 1 sensitivity limited at 
low frequencies by time-standard noise only for / ~ 1 cycle/ (10 years) or lower (Figure ([T])). 
(Within a broadband astrophysical background there may be some sources strong enough 
to be detected individually, i.e. detectable GWs coming from specific directions. Figure Q 
shows how strong an individual source would have to be for detection above the timekeeping 
noise limit. When simultaneously timing several pulsars, the directional property of the 
timing response to gravitational radiation from a single source can be used in the same way 
as for ground-based networks of broadband GWs to improve SNR, e.g. P3] and references 
therein.) Since SNR > 5 is conventionally taken as detection threshold. Figure ([T]) indicates 
that substantial SNR improvements will be required of pulsar timing array signal processing 
for reliable detections. 
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FIG. 1. Gravitational wave sensitivity expressed as strain amplitude required for SNR = 1 in a 
1 cycle per 10 years bandwidth, as a function of Fourier frequency. Sensitivity limited by various 
noises are indicated: green is station location noise [[llj; black and blue are respectively due to 
fluctuations in the zenith dry [13j and wet |14j troposphere; light blue is due to antenna mechanical 
noise for a 34m beam- waveguide station plOl I15j : dashed lines are for dispersion measure variations 
in the interstellar medium |2H [22j for 1 GHz observations after 99% and 99.99% calibration; solid 
black lines are for white timing noise with rms amplitudes of 100 ns, Ins, and Ips in a Fourier 
band ±0.5 cycles/day; dotted line is an approximate limit due to uncertainties in the masses of the 
planets [18[ [T^ ; dot-dashed line is sensitivity limited by a linear ion trap time standard j^lj • Also 
shown for reference is the low- frequency sensitivity expected for the LISA detector [25]. Red line 
and vertical bar shows the dependence and range of predicted signal strengths from an ensemble 
of supermassive black-hole binaries |29H32j 
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